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A VELOCITY-CORRECTION FORMULA FOR THE CALCULATION OF 
TRANSONIC MACH NUMBER DISTRIBUTIONS OVER 

f 

DIAMOND -SHARED AIRFOILS 
By H. Reese Ivey and KeitL C. Harder 


SUMMARY 


A velocity-correction formula is proposed for the purpose of 
calciilating, from the known Mach number distribution for a diamond- 
shaped airfoil at a stream Mach number of 1.0, Mach number distributions 
on the same airfoil at speeds from a Mach number of about 0.8 to shock- 
attachment Mach nimber. The time required to calculate these additional 
Mach number distributions is small in comparison with the time required 
by rigorous methods. The accuracy of the results for stream Mach numbers 
near 1.0 is of the. same order as the accuracy of the known Mach number 
distribution. Moreover, the results tend to become exact as the stream 
Mach number is increased toward that for shock attachment. An expression 
for the rate of change of local Mach number with stream Mach number is 
derived and an explicit equation for the drag coefficient as a function 
of stream Mach number and thickness ratio is given. 


INTRODUCTION 


The pressure distribution for a diamond- shaped airfoil at a stream 
Mach number of 1.0 has been calcxilated by Guderley and Yoshihara 
(reference l) . Calculations for a similar airfoil at four speeds between 
Mach number 1.0 and the shock-attachmient Mach number have been performed 
by Vincent! and Wagoner (reference 2). According to reference 2, 
similar calculations have been made by Cole at slightly subsonic speeds. 
These rigorous resiilts combined with reliable experimental results 
provide sufficient information for checking the acc\jracy of an approximate 
veiLocity-correction formula proposed for calculating Mach number distri- 
butions on a diamond- shaped airfoil. The concepts involved should 
facilitate the calculations for other shapes. 

The justification for the proposed velocity-correction formula is 
based upon its good agreement with existing rigorous calculations. 
Moreover, it is in accord with the general transonic similarity rule. 
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and with the special form of that rule for a stream Mach number of 1.0. 
Also, It agrees with Guderley's result that the local Mach numbers are 
constant for small departures of the stream Mach number from 1.0. 

The present paper gives the derivation of the method of calculation 
and its application for the determination of the Mach number distributions 
over a 10-percent-thick diamond-shaped airfoil in the Mach number range 
0.80 M(jo ^ 1.30. The results are compared with the earlier results 
of Vlncentl and Wagoner. A brief discussion is given of the movement 
of the shock wave along the airfoil surface as the stream Mach number 
is decreased below 1.0. 


MALfSIS MD DISCUSSION 


The velocity-correction fonniila to be presented for diamond- shaped 
airfoils is based upon the combination of two heretofore unrelated 
concepts. One concept is obtained from the special form of the transonic 
similarity irule (reference 3) when the free stream is sonic and the 
other from an examination of the transonic approximation for Prandtl- 
Meyer flow. 

Transonic similarity rule for sonic stream .- The form of the transonic 
similarity rule for a sonic stream is derived in the appendix and the 
result obtained is 


M^ - 1 = 5^/^C(x) (1) 

where M is the local Mach number, C(x) is a function depending upon 
the body shape, and 5 is a parameter (such as the airfoil thickness 
ratio or the angle of flow deflection for Prandtl -Meyer flow) used to 
differentiate bodies of the same family. The concept to be obtained 
from the similarity rule is the role played by C(x). For the purpose 
of this paper, it is sufficient to consider the flow past the nose of 
a wedge and the Prandtl -Meyer expansion around a convex corner. The 
flow deflection 5 is taken as positive for the wedge (compression) and 
negative for the Prandtl-Meyer corner (expansion). The stream is sonic. 
Equation (l) may be applied to both of these flows but C(x) will be 
different for. the two cases. For Prandtl-Meyer flow C(x) is a constant 
and for the flow past the wedge C(x) corresponds to the Mach number 
distribution given in reference 1* Thus, C(x) determines the shape of 
the Mach n umb er distribution along the surface. An illustration of 
equation (1) is presented in figure 1 to show the role of C(x). It 
should be noted (fig. 1) that all curves originate at the point M = 1, 

5=0. As 6 is varied from zero, M^ - 1 oc 6^/^ for all points on 
either surface. 



NACA TR 2527 


3 


Transonic approximation for Prandtl -Meyer flov .- The concept of 
the role played hy 5*, the amovmt the free stream must he deflected to 
reach sonic velocity, may he Illustrated hy considering Prandtl -Meyer 
flow. For Prandtl -Meyer flow. In order to reach sonic velocity the 
stream must he deflected through an angle S*, given approximately 
hy (from reference 3 "with 2 (Moo - 1) heplaced hy (M„2 - 1)): 



( 2 ) 


where M^o Is the stream Mach number and C = 
the, ratio of specific heats . The value of 5* 


j (7 + 1) 


2/3 


7 helng 


Is taken to he positive. 


Inasmuch as, for Prandtl -Meyer flow, the flow deflection may he 
considered to consist of more than one deflection (fig. 2), the ejqpresslon 
for the local Mach number may he written: 


m 2 - 1 = dh - • (3) 

where 5 Is taken as the flow deflection from the free-stream direction. 
From equation (2), 6 * = 0 when M^ = 1, and, for this case, equation (3) 

reduces to the transonic slmllsirlty rule. A plot of equation ( 3 ) Is 
shown In figure 3* It Is Important to note that the curve for M„ > 1 
Is the same as the curve for M^ = 1 except that the origin of the 
Moo > 1 curve has been shifted to the point 5*. This result shows that, 
for Prandtl-Meyer flow, the expression for m2 - 1 Is of the same form 
as the transonic similarity rule for a sonic stream, providing 5 Is 
measured from the 5* direction. 

Velocity-correction formula for dlamond-s^ped airfoil .- The 
fundamentail assumption made to develop a velocity-correction formula 
for the diamond-shaped airfoil Is that the concept of C(x) obtained 
from the slmlleirlty rule for a sonic stream and the concept of 5 * 
Illustrated hy Prandtl-Meyer flow may he combined In the form 

M^ - 1 = C(x) |b - B*|^/3 (4) 

It should he emphasized that equation (4) Is based upon an assumption 
and Is not claimed to he rigorous. Typical plots of equation (4) are 
given In figure 4. 

The C(x) used for the diamond-shaped airfoil Is obtained from the 
solution given In reference 1 for this airfoil for M^ = 1. Thus, 


Moo= 1' 


( 5 ) 
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If 1:116 stream Mach nurnher is greater than 1.0, the flow must he 
decelerated (compressed through an angle 6*) to reach sonic velocity. 
In the approximation of the transonic small -disturbance theory , the 
compression of the flow caused by a shock wave is isentropic. Thus, to 
the same order of approximation, 8* may be determined from either 
Prandtl -Meyer flow or shock-wave relations. It can be shown that 8* 
computed from shock-wave relations obeys the transonic similarity rulej 
that is. 


~ ^ _ If* 

5«2/3 ■ 


(S) 


where K*" is a constant. Note that equation (6) is very similar to 
equation (2) which is valid for Prandtl^eyer flow. The values of 8* 
used in the computations were taken from shock-wave tables to insure 
that the results obtained by the use of equation ( 4 ) would fair smoothly 
into those given by shock-expansion theory. The angle 8*, based on 
shock-wave relations, is shown as a function of Moo figure 5* 

The concept of 8* previously presented appears to have physical 
significance only for stream Mach numbers greater than 1.0. In order 
to caJ-culate flows with a stream Mach number less than 1.0, the curve 
for 8* was extrapolated by assuming that 8* was an odd function 

of Moa^ - 1. 


Conditions satisfied by proposed velocity-correction formula . - 
Equation ( 4 ) contains the transonic similarity rule as a special case. 
Eliminating 8*" from equations ( 4 ) and (6) leads to the following 
equation: 


M'" - 1 = C(x)8‘ 


1r2/3 


W 


/2 


8K*' 


2/3 


( 7 ) 

When M^ = 1, the local M^ - 1 distribution is proportional to 8^/^ 
as required by the transonic similarity rule. Moreover, when M„ and 

8 are varied so that > ^ remains constant, the expression 

for the local Mach number at the surface becomes 


M' 


2 - 1 = C(x)8^/^ 




3/2 
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This expression Is In agreement with the coirespondlng result of the 
transonic similarity rule. 


Equation ( 4 ) Is also In accord with the result of Guderley 
(reference 4 ) that = 0 for all shapes having finite thickness. 

V^VmofI 

Thus, from equation (7)^ at a fixed position x. 


and hence 


dM _ 1 

fc(x)) 


Moo^ - 1 

ClMeo " ' 

u* / 

' ir 

m2 - 1 


1/2 


(8) 



The result of Guderley that 



= 0 has heen used by Vlncentl and 


Wagoner (reference 2) to obtain the result that 


£ 24 ) 

,<3Mco/M„=l 


2 

7+1 


(fd) 


where c^l Is the drag coefficient. The method of the present paper 
gives this same slope of the drag curve at = 1. 


At a speed slightly above that for shock-wave attachment, the local 
Mach number over the front half of the diamond- shaped airfoil becomes 

sonic. Equation (8) gives =“ for this condition, whereas 

shock-wave tables Indicate that ts extremely large but finite. 

This difference between infinity and a very large quantity is of no 
practical importance In the present considerations. 

If the proposed velocity-correction formula is to yield good results 
throughout the speed range from =1 to the speed for shock attach- 
ment, the results should fair smoothly Into those for purely supersonic 
flow with attached shock waves and also into those for subsonic stream 
Mach numbers. For the case of flow with attached shock waves, 5* is 
greater than 5 and the supersonic local Mach numbers may be estimated 
by expanding around a Prandtl -Meyer corner from the 5* direction to 
the 6 direction by use of equation ( 3 ) or, preferably, the exact 
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expression for Prandtl-Meyer flow for the deflection 5* - 5. Local 
Mach number dlstrihutions were calciilated by both the present method 
and by the exact shock- expans ion method (reference 5) for a 10-percent- 
thlck diamond'-shaped airfoil in the Mach number range 1.28 ^ ^ 1.55* 

The maximum difference in the local Mach numbers predicted by the two 
methods in this speed range was less than 0.01. This agreement indicates 
that the proposed velocity-correction formula not only fairs smoothly 
into the attached- shock calculations but may also be used as a means 
for calculating such flows with good accuracy. 

The calculation for the rear half of the diamond-shaped airfoil 
at subsonic speeds is beyond the scope of this paper because of the 
presence of a shock wave on the rear surface. The location and strength 
of this shock wave seem to be strongly Influenced by viscous effects. 

The general shape of the curve of drag coefficient against Mach 
number calculated by means of equation (4) is very similar to the "possible 
interpolated" drag-coefficient curve given in reference 2 except near 
attachment Mach niimber. 


The application of the method to the front and reeir surfaces of the 
airfoil are treated separately. 

Application to front surface of diamond-shaped airfoil .- In the 
sign convention adopted, 5* is always positive ^ 1^. For a wedge, 

the flow is subsonic if 6 - 8* > 0 and supersonic if _5 - 5* < 0. 

When the local Mach numbers are subsonic, the C(x) given by equation ( 5 ) 
is used. When the local Mach numbers are siqiersonic, Prandtl-Meyer 
concepts are employed as previously discussed. 


Figure 5 shows the angle 6* for sonic flow for a wedge as a 
function of stream Mach number as calculated from shock-wave relations. 

This curve is more accurate than that given by equation (6) since K* 
is not quite constant. The difference between 6* and the local 
slope 5’ is a measure of the deviation of M^ - 1 from zero. Figure 6 
shows the variation of the M^ _ 1 distributions for the front half 
of the airfoil with the stream Mach number. These results were obtained 
by multiplying the M^ - 1 distribution for M^ = 1 (reference 1) by 

the factor I — ^ -J where 6 is now considered the wedge thickness 

ratio. Figure 7 presents the corresponding local Mach nvmiber distributions. 


Reexamination of equation (4) indicates that the proposed veloclty- 
. correction formula is equivalent to the assumption that the local Mach 
nximber distribution over the front surface of a diamond-shaped airfoil 
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is the same as that for a thinner airfoil in a sonic stream; that is^ 
5-6* may he considered as the effective airfoil thickness 6g in a 

sonic stream. Figure 8 shows, for example, that a 10-percent-thick 
wedge at a stream Mach number of I.167 has the same effective thickness 
as a 5-percent-thlck wedge in a sonic stream. A 10 -per cent -thick wedge 
at Moo = 1.278 will have the same Mach number distribution as a flat 
plate in a sonic stream; namely, M = 1. 

Application to rear surface of diamond-shaped airfoil .- For all 
Mach manbers below shock attachment, the local Mach number remains 1.0 
just ahead of the corner, and therefore the Prandtl-^eyer expansion at 
the corner remains constant. The Mach number distribution on the rear 
is considered to result from a Prandtl-Meyer expansion and the reflected 
compression waves from the sonic line. A sketch illustrating this 
Influence on the rear siirface is given as figure 9- The subsonic 
Influence (due to reflected compression waves from the sonic line) must 
decrease to zero when the local Mach number bn the front becomes sonic. 
The subsonic influence on the rear is assumed to vary with stream Mach 
mflnber in the same manner as the subsonic flow over the front and C(x) 
for the rear is proportional to 

(m 2 - Dpj, - (m2 - 

where (M^ - l)pj^ corresponds to the Mach number obtained by a Prandtl- 

Meyer expansion through an angle 25. In other words, the subsonic 

Moo = T given in reference 1 is multiplied 

to obtain 'the subsonic influence at other 

Mach numbers. The corresponding M^ - 1 and M distributions are 
shown in figures 6 and 7 , respectively. 

Figure 10 gives the pressure distributions, based upon the exact 
formula for the pressure coefficient, corresponding to the Mach number 
distributions (0.8 < M^^ < 1.3) of figure 7- Note that the pressure- 
distribution curves form a confused pattern, whereas the Mach mnnber 
distribution curves axe more uniform. This behavior indicates that 
studies of transonic flow phenomena should be interpreted in terms of 
Mach number rather than pressure coefficient. 

Figure 11 presents a comparison of the Mach number distributions 
calculated by the present method with those calculated in reference 2^for 
a 10-percent-thlck diamond- shaped airfoil. The Mac)i number distributions 


influence on the rear for 


by the factor 


2/3 
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obtained in reference 2 are presented as two curves. This choice arose 
in reference 2 because there was some question whether - 1 should 
be replaced by 2(M - 1). For the 10-percent-thick diamond-shaped 
airfoil, the Mach number is known to be exactly 1.0 Just ahead of the 
corner and 1.1(85 Just after the Prandtl-Meyer expansion at the corner 
for inviscid flow. 

Drag coefficient .- Figure 12 shows the variation of drag coefficient 
with stream Mach number for the 10-percent-thlck diamond-shaped airfoil. 
Separate ciirves are presented for the front and rear surfaces in order 
to emphasize the change in relative importance of the drag contributed 
by the two parts as the Mach number varies. An attempt has been made to 
reproduce faithfxilly the peculiarities and breaks in the curve, especially 
as the local Mach number on the front surface becomes sonic. The 
present results fair smoothly into the exact supersonic results; in 
fact, the supersonic results can be obtained to a high degree of accuracy 
by the present method. 

Figure I3 presents a comparison of the drag curve given by the 
proposed method with the "possible interpolated curve" of reference 2 
for a 7 ‘87-percent- thick diamond- shaped airfoil. The "boxes" shown 
in figiire I3 represent the choice of drag coefficients which are due to 
the choice of press\ire coefficients and stream Mach numbers presented 
in reference 2. The two cvcrreB are very similar except near attachment 
Mach number. Reference 2 shows a fairly small slope for the drag curve 
in this region, whereas shock-expansion theory and the present method 
Indicate that this slope becomes extremely large. 

By use of the present concepts an explicit expression may be derived 
for the drag coefficient of a diamond-shaped airfoil for M;„ ^ 1. The 
drag coefficient may be expressed in terms of the airfoil thickness 
ratio 6 and average pressure coefficient P as follows: 

Cd = 

Except in the equations to follow, the exact pressture relations have 
been used throughout this paper. If the exact expression for the pressure 
coefficient were used in the analysis to follow, the formula obtained 
for the drag coefficient would yield the drag curve shown in figure 12. 
However, for the sake of sinplicity, the following approximate relation 
obtained from Bernoulli’s equation and limited to Mach numbers in the 
neighborhood of unity is \ised: 


(m2 - 1) - (M<2 - 1) 


2 

7 + 1 


P 
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Then 


2B 

y-n 


(m2 - - C(x)r( 6 - 5*)^/^ + C(x)p(5 - 


— r ^ 

'3 ■' 
/_ • 


where C(x)j. = for the front stjrface at Moo = 1 and 


62/3 


C(x)jj = 


(m2 - 1)^ - (M^ - 1)^ 

r2/3 


for the rear surface at Moo = 1. 


For Prandtl-Meyer flow, M^ - l = ^3 6 ^ and replacing 

C(x)p and C(x)j^ in terms of the drag coefficient for M^ = 1 leads 
to the following approximate expression for the drag coefficient: 


=d = 


28 


5/3 


(r + 1)^/3 


1 2/3 + 


28 


'do 


,5/332/3 


(r + 1)^/3 


1 1 ““ 

” bV K*-. 


- 1 


, 3 /^ 


2/3 


( 9 ) 


where c^^ is the drag coefficient for M^^ = 1. 


Equation (9) is not sufficiently accurate for general use because ' 
of the rather severe limitations of the approximate pressure formula 
hut is xiseful for illxxstratlng the variation of drag coefficient with ? ■- 
stream Mach number. 


Equation (9) separates the dreig coefficient for a diamohd-shaped 
airfoil into two terms; (l) a supersonic term contributed by the Pran^*l 
Meyer corner which remains constant and (2) a subsonic teimi which 
decreases to zero in the vicinity of shock-attachment Mach number. 

At stream Mach numbers less than 1.0, the drag qonti4buted by the 
supersonic term is no longer constant but is decreased, because a shook 
wave moves forward along the rear surface. This forVard' movement 
decreases the extent of the supersonic region which /is gradually replaced 
by the symmetrical subsonic pressure distribution. ,The existence of a 
shock wave on the rear surface at these speeds makes the calculation of 
the pressure distribution very difficult since the ;iocation of the shock 
wave is strongly affected by boundary -layer , separation. For this reason, 
the total drag coefficients in the subsonic range are not presented. 
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CONCHTDIBG REMARKS 


A veloci'ty-correc'fcion formula has "been proposed for "the purpose of 
calculating, from the known Mach numOber distrihutlon for a diamond- 
shaped airfoil at a stream Mach number of 1.0, Mach nximher distrihutions 
on the same airfoil at speeds from a Mach number of ahout 0.8 to shock- 
attachment Mach number. The formula exhibits the following properties: 

( 1 ) The formula contains the general transonic similarity rule as 
well as the special form for a sonic stream. 

( 2 ) The local Mach number over the front s\arface of the airfoil 
is 1.0 at the correct stream Mach number. 

(3) The drag coefficient is that given by Guderley and Yoshihara 
when M<a = 1 . 

(4) The rate of change of drag coefficient with stream Mach number 
at Mj„= 1 is the same as that given by Vincenti and Wagoner. 

( 5 ) Calculated Mach number distributions are in agreement with 
calculations based on shock-expansion theory above shock-wave attachment 
Mach numbers. 

( 6 ) The formula reduces to the proper form for Prandtl -Meyer flow. 

Pressure and Mach number distributions for a 10 -percent -thick 
diamond-shaped airfoil are presented for the Mach number range 
0.8 < M„ < 1 . 3 . The variation of local Mach number distribution with 
stream Mach number is regular, whereas the pressure-distribution curves 
form a confused pattern. This behavior indicates that studies of 
transonic flow phenomena should be Interpreted in terms of Mach number 
rather than pressure coefficient. 

By the use of the velocity-correction formula proposed in the present 
paper, an approximate expression is derived for the drag coefficient 
of a diamond- shaped airfoil as a function of stream Mach number and 
thickness ratio. The drag coefficient is separated into two terms: (l) a 

supersonic term contributed by the Prandtl -Meyer .corner and (2) a subsonic 
term which decreases to zero in the vicinity of attachment Mach nuniber. 

In its present form, the proposed velocity-correction formula 
appears to be applicable only to diamond-shaped airfoils at transonic 
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speeds and will probably require modification if it is to be applied 
to cxirved airfoils. Presumably, the present concepts may be used to 
obtain coiresponding relations for a cone. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va., August I5, 1951 
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APEEKDIX 

DERIVATION OF TRANSONIC SIMILARITY RULE FOR A SONIC STREAM 

Ttie resul'ts presen.'ted in reference 3 Lave shown ■that, when a 
perturhation velocity potential 0 is defined hy 

u 

V 

wtiere a* is the velocity of sound for M = 1, and x and y are 
Cartesian coordinates the transonic approximation to the differential 
equation for the flow of a con^iressible fluid may he written as 

^ ^xx ~ ^yy - ^ 
a* 

In addition, when 2(M - l) is replaced hy M^ - 

= («) 


= a* + 0J 

= 0y 


(Al) 


Also, from reference 3, the appropriate boundary conditions for 
equation (2) when the stream is sonic are 

0y(x,O) ^ 



m 


and at infinity 


= 0 (A5) 

where 8 denotes the airfoil thickness ratio, c is the chord, and h 
is a function describing the body shape. 

A solution 0i(xi, yi) of equations (A2), (a 4), and (A5) is assumed 

to be known. Two flows are considered to be similar if a solution 
02(^2^ Yg) satisfying equations similar to (A2) , (a 4) , and (A5) can be 
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related to In particular, similarity will exist providing A, B, C 

can be so determined that 


02 = -^^^1 
X2 = Bxj 
2^2 ~ 

if, in addition, and ^2 satisfy equation (A2) , and 

0) = 

0) = 52h(^'j 


(A6) 


(A7) 


and at infinity 




With full generality, c^' may be taken equal to C 2 since 

both bodies are in the flow field which extends infinitely far in every 
direction. If c^ = then B = 1. Flows involving different values 

of 7 need not be considered in this analysis, which is primarily 
concerned with the result for a particular gas - namely, air. 

If both 02. 02 satisfy equation (A2), it is easily 

found that 

AC^ = 1 (A8) 


From equations (a 6 ) and (AT) 

= v(|) = £ 02y^(x,O) = £ 52 h(|) 

from which 

2. = 

A 82 


(A9) 
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From eg^tions (a 8 ) and (A9) , 


'■(S) 


1/3 


A = 


_ 


(AlO) 

(All) 


From equation (A3) 

!< 2 ^(==, 0 ) = - 1 ) = A0i^(x,O) = A 5 ^ (Mi^ . ( 102 ) 

and, from equations (All) and (A12), 




2/3 Ml - 1 


Inasmuch as 0]^ vas assumed to he the known solution, the value 


of 


- 1 


8 


2/3 


at the surface may. he replaced hy C(x), a known function. 


and the transonic similarity rule for a sonic stream yields the result 
that, on the surface, 


- 1 = 6 ^/^ C(x) 


(AI3) 
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Figure 3-- Prandtl-Meyer flow plotted to illustrate role of 


- 1 = Cl5 - 6*1 2/3; 



6 *. 
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Figure k.- Illustration of proposed velocity-correction formula showing 
combination of concepts of C(x) and B*. - 1 = C(x)]B - 5*1^/^; 

- l\3/2 







Figure 5 .- Variation of 6 * vith stream Macti nmber 
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Figure 6 .- Calc-ulation of - 1 distriimtlon for 10 -percent -thick 

diamond-shaped airfoil. 
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Figure 7 -- Local Mach, number distrlbutiona on lO-percent-thlck diamond- 

shaped airfoil. 
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Figure 10.- Pressure distritiutions on a 10 -percent-thick diamond- 
shaped airfoil. 
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rigure 11.- CoBUparlaon of Macli nuuiber dlBtrl'butioiis on a 10 -percent -thick 
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